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ABSTRACT 


This paper considers the irrotational motion of an incompressible, inviscid fluid con- 
tained in a partially filled tank. The tank is subjected to both transverse and rotational 
vibrations whose frequencies are near the first natural frequency of small free-surface 

O 

oscillations. Following a method suggested by Hutton , the analysis is performed re- 
taining higher order terms in the free-surface dynamic and kinematic boundary condi- 
tions. The theoretical investigation predicts the forcing frequency ranges, for various 
combinations of rotational and translational motion, over which there are stable, 
steady-state, harmonic, planar and nonplanar motions. The least stable case occurs 
when a combination of motions occurs in mutually orthogonal planes. 
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SECTION 1 


INTRODUCTION 


The linear theory of the small oscillations of a free surface in a gravitational field 
appears in essentially complete form in Lamb's Hydrodynamics-^. The advent of the 
missile age has stimulated renewed interest in this problem, and in recent years 
numerous papers have been written on specific problems in this area. 

The dynamic response of the liquid propellant in the tanks of a space vehicle can affect 
the stability of the vehicle; this influence can be alleviated in many ways, among which 
are the proper choice of tank form, of tank location, or the introduction of baffles into 
the tank. 

These fluid oscillations, resulting from such sources as perturbation of the trajectory, 
have been shown experimentally to be most critical when the excitation frequency is in 
the region of a natural frequency of lower mode fluid oscillations. 

Eulitz and Glaser 2 have compared experimental results with the previously obtained 
theoretical solutions, which are obtained from a linear boundary value problem. Within 
the framework of linear theory, the free surface of the fluid in a container undergoing 
transverse harmonic vibrations should exhibit a steady-state, planar, harmonic motion 
at all frequencies except resonance. Eulitz and Glaser claimed thorough agreement 
between the experimental results and the linearized theory. 

Hutton 3 notes that the free surface of a fluid in a container undergoing transverse har- 
monic vibrations does not necessarily exhibit a steady-state harmonic motion. In fact, 
if the container is excited at a frequency well below the lowest natural frequency, p^j, 
of small, free-surface oscillations, the steady-state fluid motion is harmonic with a 
constant peak wave height and a single nodal diameter perpendicular to the direction of 
excitation. The wave height increases with an increase of the excitation frequency. 
When the excitation frequency is close to but smaller than p^, the smoothly oscillating 
free surface changes to a violently splashing condition. As the frequency increases, 
this motion continues until a frequency greater than p^ is attained. Additional in- 
creases in the excitation frequency reduce the wave height up to the point where the 
cycle begins again as the next resonant frequency is approached. 

Hutton shows that the sloshing motion can be accurately predicted in an inviscid liquid 
if the analysis includes appropriate nonlinear effects. 

The present paper is an extension and generalization of Hutton’s work to include not 
only transverse harmonic oscillations of the container but also rotational harmonic 
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oscillations. Again, the appropriate nonlinear effects are included. A comparison is 
made between the two solutions and, in particular, the stability of the nonlinear motion 
is studied. 
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SECTION 2 


DEFINITION OF THE BOUNDARY VALUE PROBLEM 


The problem under consideration is that of a tank, partially filled with a nonviscous, 
incompressible liquid, which is mounted in a system which is moving along a prescribed 
path. Perturbations of the path of the system cause the liquid to oscillate. There exist 
two possible types of motion that should be considered. The first is that of surface waves 
of large amplitudes, possibly of low frequency, which could actually damage the tank 
structure. For the most part this type can be controlled by suitable baffles in the tank. 
The second type, which will be considered here, is that of surface waves of small am- 
plitudes with a frequency near the natural frequency of the control system on the tank, 
i.e. , the natural frequency of the liquid-tank configuration. 

Since the tank is in motion along some path, it seems reasonable to refer its motion to 
an inertial coordinate system, for example the earth. However, if any type of measur- 
ing device is attached to the tank, then it measures quantities in terms of a tank-fixed 
reference frame which is moving relative to the inertial system. Thus, it is necessary 
to be able to express the tank -fixed system in terms of the inertial system and vice 
versa. 

Let Yj be an inertial Cartesian coordinate system with origin O' and coordinates y^; 
and let be a Cartesian coordinate system moving relative to Yj, with origin O and 
coordinates Xj. Then, instantaneously, the position of a particle moving with the Xj 
system can be described in the Y^ system by 

y. = Z. + a., x (2.1) 

i i Ji 1 

where the summation convention is being used and Latin subscripts take on the values 
1, 2, and 3. In (2.1), Z.(t), with components measured in Y^ give the instantaneous 
displacement of O relative to O'; and 

a ijO^ = cos(x., y.) (2.2) 

measures the instantaneous rotation of x^ with respect toyj. Subsequently the follow- 
ing notation will be used: a barred vector has components measui’ed in Yj and an 
unbarred vector has components measured in X^. 


Since ajj are a set of direction cosines, they satisfy, for any t. 


a 


ik a jk 


6 .. 


(2.3) 
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where 5^ is the Kronecker delta. Denote df/dt by f and take the derivative of (2. 3); 
then 


a a + a a =0 
ik jk ik jk 


(2.4) 


Define 


oo.. = a., a., 
ij ik jk 


(2.5) 


Thus , using (2 . 4) and (2.5), one has 


co.. = a a 

ji jk ik 


- a a 
jk ik 


■ a a — - oo.. 
ik jk ij 


( 2 . 6 ) 


That is, (Ojj is a skew- symmetric second order quantity which can be shown to be a 
second-order tensor. 

A dual vector oOj can then be defined such that 

oo.. = -e... oo, 

i] i]k k 

where is the third-order alternating tensor. Thus, from (2.6), 


a ik E jk ^ijk^k 


(2.7) 


where 60^ is the angular velocity of Xj with respect to Yj measured in Xj. 

The absolute velocity of a particle whose position is described by (2. 1) can be found by 
differentiating (2. 1) with respect to time. This operation gives 


y = Z. + a x. + a..x = q. 
i i ji J l 1 j i 


( 2 . 8 ) 


where q^ is the velocity measured in the Yj system. The measuring device fixed on 
the tank measures q^, where qj is the velocity measured in the Xj system and 


<L = a - q i 
1 ij j 


(2.9) 


Using (2.8), q^ can be written as 

q i * a ij( 2 j + “kj^k + \j\) 


( 2 . 10 ) 
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V 


With (2.3), (2.7), and the fact that Zj = ay Zj, (2.10) becomes 
q. = Z. + x. - . 00 . x, 

l l l ikj j lc 


or 


q. = Z. + x. + e... to. x. 

l l l ijk ] k 

where the skew- symmetric property of the alternating tensor has been used. 
The next quantity of interest is the absolute acceleration 


a. = q. 
x 1 


From (2.9), (2.3) can be written as 


d 

dt 


d q ; 


(a. . q. ) = a.. —1 + a. . q. 
\ ji ]/ ji dt ji j 


The quantity of interest is a^ = a- aj; thus from (2. 14), 


dq. 

a , = a , - a -- + a . - a --q. 

k ki ji dt ki ji j 


or 




\ dt ^ijk^i^j 

Here it is noted that the velocity q^ is a function of not only the time, but also the 
dinates x^, which are also functions of time; thus 

dq k aq k 3q k 

+ 


dt at a x. dt 

i 


From (2. 12), 
dx. 

l 

dt 


x. 


. = q. - z. - e. co.x, 

l l ijk j k 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


(2.14) 


(2.15) 


(2.16) 

coor- 


(2.17) 


(2.18) 
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Then (2.18), (2.17), and (2.16) give, finally, 


a = + ^ 
k B t S x 


-fq. - Z. - e.. co.x ] + e... co. q. 
I i i ijs j sj ijk i j 


(2.19) 


In vector form, (2. 11) and (2.19) appear as 


q = q + co><r + r 
o 


( 2 . 20 ) 


and 


|y+o)Xq+^q-q o -(oXrj-V q 


( 2 . 21 ) 


where q Q is the velocity of O relative to O' 

The Eulerian equations of motion for an incompressible, inviscid fluid are, in the 
inertial system, 


a. = F. - 1|P- 

i i pay. 


( 2 . 22 ) 


where F. is the specific body force, p is the density, and p is the pressure. Since 


y. * y,fx„ x , x , 
J i\ 1 2 3/ 


3 p Bp 


3 \ 


= a. 


BP 


B y 9 \ 9 y i ki B x 


(2.23) 


Transform (2.22) to the tank fixed system using (2.23); then 


_ 1 Bp 

a = F - — a . a, . 

j j p ji ki^ 


(2.24) 


or 


a = F - 

j j pax. 

Assuming that the motion is irrotational , there exists a potential, $, such that 

30 


q. 


Byj 


(2.25) 


(2.26) 
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II 


The incompressibility assumption implies that 


Equations (2.26) and (2.27) then lead to 



Transform the above to the tank fixed system, noting that 


(2.27) 


(2.28) 


9 0 __ 9 0 ^ _ 9 0 

" a,dJ \ 

Then 

9 0 

q = 

9 x. 

] 


9 x. 
J 



(2.29) 


(2.30) 


(2.31) 


Thus, the solution of Laplace's equation furnishes a possible potential function for an 
incompressible, irrotational flow. In order to determine exactly which potential func- 
tion is the solution, certain boundary conditions need to be prescribed. 


Consider a tank of arbitrary shape partially filled with fluid. Assume a constant accel- 
eration is acting along the X3 axis. The surface of the liquid then assumes a planar 
surface normal to this axis, this surface being called the free surface or quiescent 
free surface. The origin of the Xj system is taken at the center of gravity of the accel- 
erating fluid system. The motion of the tank-fixed system Xj relative to Yj, charac- 
terized by Zj and 00^, are oscillatory motions superimposed on the constant-acceleration 
motion. These motions will induce perturbations or disturbances of the free surface. 
The measuring device traveling with the tank sees only the forcing motions or 
perturbations . 
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In this analysis it will be assumed that the tank is rigid. With this in mind, the bound- 
ary condition on the wetted surface of the tank must be that the velocity of the liquid 
normal to the tank wall must equal the normal component of velocity of the tank itself. 

Thus, if tv is the unit exterior normal to the tank and qi = - , 

<3 xj 


-•'iH; = ^i[ z i + t ijk‘°i !t kl (2 - 32) 

where = 0 for a rigid tank. 

There are two conditions at the free surface. Denoting the disturbed free surface by 
77 (Xj , x 2> t) and the unit normal to the quiescent free surface by n i; the kinematic con- 
dition that a particle of fluid which travels with the free surface as it moves must have 
the same velocity as the free surface itself is given as 



(2.33) 


where x 3 is the displacement of a particle in the x 3 direction, and, as in (2. 17), 


d 

dt 


_a_ 

at 


+ 


^kj_ 

at ax, 


(2.34) 


Expanding (2.33) and using (2.34) and (2.12) 


St? a-n . ^ a r? . „ 

— - + — - x, + — - x„ = q - Z„ - w.x, , onx„ = 7 ) 

at ax i ax 2 M 3 3 3 jk j k 3 ' 

-L A 

But since n^ = n 2 = 0 and n 3 = 1, the right hand side of (2. 35) can be written as 

(s ‘ z i - f ijk“j\) n i 
Thus (2.35) becomes 


(2.35) 


+ i. - - i. - - u.x, )„ , 

at ax. 1 \ 1 1 13k ] k/ i 


on x g =77 


(2.36) 


The second condition at the free surface is a dynamic one which states that the pressure 
at the free surface of the fluid must equal the ambient pressure. To find the form of 
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this boundary condition it is necessary to integrate the equations of motion. Substitute 
(2.19) into (2.25); the equations of motion become 


at 


a x. 

i 


q. - Z. - e 00 x + e co. q. 
L x i its t sj ljk i j 


F k- 


i ap 
p sx k 


(2.37) 


If the only specific body force is that due to the gravitational field in which the tank sys- 
tem is moving, (2.37) can be integrated directly: 


P - P, 


= -[ + - Zi)Mv“j X k) q l 'If 


(2.38) 


where p Q is the ambient pressure and a is the magnitude of the acceleration of the tank 
system. It is assumed here that p Q is a constant. Thus, the second boundary condition 
at the free surface is 


a 0 

aT = 


1 /JL 0 

2 \3x. 



+ Wj x k-|^-,on x 

ljk J a Xi 2 


V 


(2.39) 


In summary, the mathematical description of the motion of an incompressible, irrota- 
tional fluid confined in a moving, partially filled tank subject to translational and 
rotational perturbations is 


2 

V 0 


0 


(2.40) 


a 0 
v i a x. 

i 






(2.41) 


on the wetted surface, where y ■ is the unit exterior normal to the tank; 


12 + 12 Ik 
at ax. i 

X 


/ d 0 

\ax. 


+ Z. + e. oj 
i i]k 


A)"!' 


on x_ 


V 


(2.42) 


where n^ is the normal to the free surface; and 


a 0 

a~t 


at} 


1 1 a 0 
2 lax. 



90 

f ijk aJ j \ ax. 


on x 

3 


V 


(2.43) 


If the free surface oscillations are sufficiently small, terms of second order in the 
velocities can be neglected. Then (2.40) to (2.43) become 

y 2 0 = 0 (2.1a) 
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-U = v. [z. + e. (o.x. 1, on S 
1 3 x. ill i]k ] k J o 


(2.2a) 


where is the wetted surface. 


dr) = / 30 

b t lax. 


(■—■ + z - + c... to. X^)n. 
VBXj i i]k ] k/ l 


, onx„ = 


(2.3a) 


3 0 

— = an, onx 3 = 7) 


(2.4a) 


The last two equations can be combined into a single condition: 


2 

13 0 ( 30 ' \ 

— = - ( — — + Z. +€■.., oj. x, n., on x„ = r) 

a 2 \ 3 x. i i]k 3 k/ l 3 

o t 1 


(2. 5a) 


If is is assumed that Zj and coj can furthermore be represented as harmonic oscillations, 


<°>eV*t 


Z. = Z) 'e 
l l 


(0) wt 

o. = CO. 


(2.6a) 


then it may be assumed that 


0 ( x i’ l ) = ^i)®^ 


(2.7a) 


The problem then reduces to 


V ip = 0 


(2. 8a) 


on S 

i 3 x. ill ljk j k j 6 


(2.9a) 


+ z (°) + t = „ 

a \Bx. l ijk j k/ i 3 


(2.10a) 


The solution to (2. 8a) through (2. 10a) can always be obtained if the tank is a prismatic 
cylinder with x 3 parallel to a generator and with cross section such that ip is separa- 
ble in the appropriate three-dimensional coordinate system. 
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SECTION 3 


NONLINEAR SLOSHING IN A CIRCULAR CYLINDRICAL TANK 


The natural course here is to express the problem in terms of cylindrical polar coor- 
dinates (r, 0, z). Before doing this, consider (2.41); there exist two segments of 
wetted surface: the side of the tank and the bottom of the tank. Considering the side 
first, the normal pj has the following components: 

u 1 = cos 0 1 


v 2 = sin 0 


v 


3 


0 


Thus, (2.41) becomes 


(3.1) 


-u .. v = v Z + v Z + 1/ e mx + v <r co x. 

lgx^ 2Sx 9 11 2 2 1 ljk j R 2 2jk j k 


(3.2) 


On the bottom of the tank the normal v ^ has the components 


^ = ° 


, 2 = 0 


V = -1 
3 


(3.3) 


Thus, (2.41) becomes 


B 0 


B x 


3 



C 3 j k t i\ 


(3.4) 


Note that Z 3 is absent in (3.2) and that 0O3 is absent in (3.4). This condition will in fact 
be the case for any cylindrical tank whose generators are parallel to the x 3 axis and is 
not merely a peculiarity of the circular-cylindrical tank. 


In the following, (2.40) through (2.43) will be transformed into cylindrical polar coor- 
dinates. So far all the quantities in these equations have been measured with respect 
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to the tank-fixed rectangular Cartesian coordinate system. These quantities may also 
be expressed in a tank-fixed, cylindrical, polar coordinate system. In doing so, it is 
also convenient to shift the origin from the center of gravity of the fluid to the geometric 
center of the quiescent free surface. Let 


X 2 = y 


x = z 
3 


and use the usual transformation equations to cylindrical polar coordinates 


x = r cos 0 


y = r sin 0 


(3.5) 


Then (2.40) becomes 


2 

V 0 


2 2 2 
8 0 + 1 8 0 + 1 3 0 + 8 0 

2 r 8 r 2 2 2 

8 r r 8 0 8 z 


= 0 


(3.6) 


Let 


b. = e... co. : 

l - - 


"ijk t °j\ 


b = 00 z - go r sin 0 
1 3 O 


b 2 = °°3 r C ° S 9 ~ z 


b 3 = oo 1 r sin 0 - r cos 0 


(3.7) 


Denote the quantity Z| by u^; (3.2) and (3.4) become, respectively, 


— = - u , cos 0 - u sin 0 - z go cos 0 + z go. sin 0 
8 r 1 2 2 1 


(3.8) 


d 0 . . 

- — = - u„ - co, r sin 0 + go„ r cos 0 
dz 3 1 3 
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(3.9) 



I 


It may be noted in (3.8) that 003 is absent. This is only true for a circular 
cylinder. To see this, consider the last term, tOj in (2. 41). In vector form 

this is 

V . (00 X p) = co . (px v ) (3.10) 


where 


p = ( r cos 0, r sin 0, z) . 

On the side of the circular-cylindrical tank, Ui has components given by (3.1). Thus, 
the third component of p x p is 

(px V) = r cos 0 sin 0 - r cos 0 sin 0 = 0 (3. 11) 

O 

Thus, from (3. 10), 0)3 is absent in the expression i^ieijkCOjXj^ for a circular cylindrical 
tank. Since this term in (2. 9a) is the only place enters, 0)3 is absent from this 
boundary condition for this special tank configuration. 


On x 3 = 17 (r, 0, t), the unit exterior normal has the following components: 



(3. 12) 



Therefore, from (2.11), (2.42), (3.7), and (3.12), one of the boundary conditions at 
z = t) becomes 


9 0 
9 z 


- u 3 - r sin 0 + 


co 2 r cos 0 = 


9r) 

FT 


9t) 9 0 

FF Fr 


1_ 9tj 90 
2 90 90 " U 1 


(ose fa 


sin 0 


12)- u (sine^a 
ee/ 2 \ ar 


cos 0 


977 \ 
90 / 


+ U 3 -.T] (si 


sin 0 


9 77 | cos 0 
9 r r 


in 

90 


)-<o 2 „(cose fa 


sin 0 9r)\ 9rj 

r Fez “ Fs Fe 


(3.13) 


The other boundary condition at z = 77 can be obtained from (2.43) using (3. 7) and (3. 12): 


90 _ 1 

Ft - + 2 


/90 

) 2 +(1MV 

\9r 

/ \ r 9 0 / 


290 
9 r 


^ cos 0 + u 2 sin 0j 


290/ „ \ „ 90 2 2 2 

( u cos 0 - u sm 01+2 — u„ + u + u + u 

r 9 0 V 2 1 / 9 z 3 1 2 3 


3-3 



(3.14) 


3 0 

3r 

/90 . 

'^W 81 


03 


sin0 + 


2 30 , 

~ su 

0 

cos 0 3 0 \ 
" r 3?/ 


o® _ & o© . . 30 90. - „ 

+ “2 a 7 ” cos 9 ' T ” sm 6 + “s + r aT ( “i sm 9 ‘ “2 cos 9) 


Thus, if the cylindrical tank has radius a, and the original depth of the fluid is given by 
z = -h, the motion of the contained fluid in a circular cylindrical tank is governed by 
the following partial differential equation, subject to the given boundary conditions: 


for 


2 2 2 
3 0 _1 3 0 _1_ _3 0_ + 3 0 

2 r 3 r 2 2 2 “ 

3 r r B0 3z 


0 ^ r < a 
0 < 0 £ 277 


and 


On r 


-h < z < T] 


= - u cos 0 - u sin 0 - zw cos 0 + z co, sin 0 
3 r 1 2 2 1 


On z = -h 

3 0 


— = -u„ - go, r sin 0 + co n r cos 0 
c) z 3 1 2 


On z = rj 


30 . A ^ Q 37? 377 30 1 3?7 90 

— - - u - go, r sin 0 + go r cos 0 = — - - — — - ~z tt 

az 3 1 2 3 1 3 r 3 r 23030 


- u. 


- u. 


cos 


sm 


. 377 sin0 3t?1 I . . 3 7? cos 0 3 7? I 

3r r 3 0J l'l 3r r 3 0 J 

3r? cos 0 3 7 ?] [ „ 37? _ sin0 3 t? 

9 37 + r 30J CO 2 TJ l COS0 3r r 30 


3 7 ? 
! 3 3 0 


30 , 1 

at = av + 2 


(£*\ 2 + (i m ) 2 + (^V 

\3 r / \ r 3 0/ \3 z / 


(3.15) 


(3.16) 


(3.17) 


(3.18) 


+ — u, cos0+u„sin0 
3 r 1 2 


1 30 f . . .1 A 30 If 2 2 2] 

u cos 0 - u sin 0 + — — u„ + — u, + u n + u I 

r 30 12 1 J 3 z 3 2 [ 1 2 3 j 
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[ 30 „ 1 30 . 1 30 30 _ . „ 

— COS0- - — smej + to 3 — + t^- r [tOj sm 6 - u> 2 cos 0] 

f 3 0 


. . cos 6 3 01 
sin 9 + — 

r 30 J 

To simplify the problem described, make the following transformation: 

0 (r, 0 , z, t) = 0 (r, 0 , z, t) + u r cos 0 + u r sin 0 - z r m sin 0 

+ Z OJg r cos 0 

Equations (3. 15) through (3. 19) become: 

v 2 0 = 0 

On r = a 

3 0 


(3.19) 


3r 

On z = - h 


= 0 


30 

= - u - 2 r O) sin 0 + 2 r« cos 0 
3z 3 1 2 


On z = r) 


3 0 


- u i r cos 0 - u 2 r sin 0 + r] rco^ sin 0 - t)cu 2 r cos 0 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


and 


aT1 2 


a an) ♦&) aaan 


3 2 / 2 2 2 2 \ 30 / u 3 , 

+ 2 r sin 0 + co 2 cos 0j + 2 ^-^rco^ sin0 - rco 2 cos0 + — J 

2 

- 3 r oj i o) 2 sin0 cos0+ r sin0 - u 3 co 2 r cos0 - 770 J 2 u 1 - u 2 o) 3 r cos0 
+ 7)0^ + 77co 3 co 2 rsin0 + u^rsin© + w 3 + r cos0 (3.24) 


- u „ + 2 rm cos 0 - 2 ru> sin 0 = | U - - _L w (3.25) 

3 z 3 2 1 3t 3r 3r 2 30 30 w 3 30 


3-5 



The free surface rj is one of the unknowns, but it may be eliminated between (3.24) and 
(3.25) as shown in Appendix A. Equations (3.24) and 3.25) may thus be replaced by 
(A. 12): 


where 


a b 
11 00 , 

a + a - + a 

00 01 a 


02 a 


ab ab abb ab 
11 01 12 00 11 11 00 22 00 
+ + + 


a 


O' 


2 a 


+ B 11 +B 12 +B 13 " + o:” 2 [rdj u 2 r sin 0 


- b u r sin 
11 2 


2 2 2 
r cos 0 


„ . , 7 -l|3 22.2 3 

0 + rcr 1 b Q0 ] + « % sln 9 + 2 

+ 2 £ r r^ u g - 0' 1 J - 3r 2 oo i uJ 2 sin 0 cos 0 + 2u 3 ro 1 + U 1 W 3 r sin 0 


+ “ 3^0 


- bJ*- 1 

[b 

/ 

22 1 

l 00 

-If, 


3 2 2 

. 2 

+ a b 

01 

3 

u 

1 

+ 

sm 

2 1 



J - a _2 [b n - raj[l 


00 2 


- 3 r co 1 a> 2 sin 0 cos 0 + 2u g rc^ + r sin 0 + u> 3 £ 


-1 . 

*12“ V 


-1 

' 23 ^ "00 ~2 


22 -2 / 

+ a,„a u„ r sin 0 + B„ 0 a *|b_ + u„ r sin 0 | + — a r sin0 


.22 2 
+ u r sin 
2 


T} 

2 \ 33 -2 r 1 2 / 4\ 

in e) + —a [b 00 + ursine] + 0(77 J =0 (A. 12 ) 


«8 = H (r ’ 01 °' t) 


■ !!< r ’ 9 ' »• *> 


i = |^(r ,0 , 0, t) 
^r 9 r 
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ay, by, and Bjj are functions of the potential ip and its partial derivatives, all evaluated 
at z = 0 (see Appendix A); and 

= u> 2 cos 0 - co 1 sin 8 


This replacement leads to the boundary value problem consisting of (3.21), (3.22), (3.23), 
and (A. 12), which involves only the potential function ip and the prescribed tank dis- 
placements. The tank displacements are assumed to be 


x i (t) = Cq sin to t 
0. (t) = 0^ sin cot 


-v 


(3.26) 


i = 1, 2, 3 J 

with e* and 0* "small" and co close to or equal to the lowest natural frequency p-^ 
given by 


P 


11 



tanh \ h 


(3. 27) 


where A is the first non-zero root of 
J i( A ln a ) ■ ° 

Here the x^ (t) correspond to translational motion, and the 0j (t) correspond to rotational 
motions. Since and 0* are small, it can be effectively assumed each set is the same 
for all i, say and 0 q, respectively; and furthermore it can be assumed that 



The tank velocities x i (t) and 8 i (t) are 
x. (t) = e cos co t 

f (3.29) 

0 . (t) = — cos CO t 


where 


e = coe Q 
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Following Hutton 3 , a steady-state harmonic solution to this boundary value problem is 
posed in a perturbation form, in analogy with the Duffing problem^’ ® in terms of 
the parameter e 


>P 


_1 

= [$2 (r, t) cos co t + (r, t) sinto tj 

2 _ 

+ + ip 2 ( r ) cos 2 604 + X 2 ^ sin 2tot ] 


+ e (r) cos 3 oo t + (r) sin 3 u) t 


(3.30) 


where the functions di and X for each value of n, each satisfy 
n n 


v 2 $ = 0 


3<i> „ 

= 0, on r _ a 

3 r 


3 $ 

= 0, on z = - h 

3 z 


> 


J 


(3.31) 


Here r means dependence upon r, 0, and z. A set of normal modes of vibration which 
satisfies (3.31) identically is 


A (t) cosm© 
mn 


+ B (t) sin m 0 
mn 


J (x r) 
mV mn / 


coshX (z + h) 
mn 

coshX h 
mn 


(3.32) 


where the J m are Bessel functions of the first kind of order m, for m a positive integer 
or zero; and X mn are an infinite set of numbers for each m obtained from the equation 

J '(x a) = 0 (3.33) 

mV mn / 

The functions A mn (t) and B mn (t) will be called the generalized coordinates of the mn'th 
mode; they depend only on the time, t. The natural frequency of small, free-surface 
oscillations in the mn'th mode is denoted by P mn - When the tank displacements are 
harmonic motions at a frequency close to or at the lowest natural frequency, p-Q, 
associated with the mode, the generalized coordinates A^ and B-q dominate all 
other generalized coordinates. Thus, it is assumed that the first order terms, ipi and 
Xj, in (3. 9) contain only the Ji mode; thus ipi and Xi are chosen as 
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cosh 


i , , N - ’(\i< z + h > 

*1 " |t 1 < T)oos9 + t 3 (T)S “ 9 l J 1 ( X 1 l r ) coshx h ' 


(3. 34) 


cosh 


X 1 = 4 f 4 fr)si»9]j 1 (x il r) 


(3.35) 


where the transformations 


1 2 

= ~ € 3 to t 


11 


2 , 

CO 1 


- ve 


(3. 36) 


have been used. In (3.36), v is a dimensionless measure of frequency and r is a dimen- 
sionless time parameter. A derivation of the above transformation is given in the 
appendix. 

As shown in the appendix, (3.30) is then substituted into (A. 12). Equate the coefficient 
of <r^to zero: 


(n>P lz ~ P ii^i) COS u)t + ( aX 2z ~ l i 1 | X 1 ) sinwt = 0, on z = 0 (A. 22) 

Thus (A. 22) is satisfied identically for all time if j/q and Xj are chosen as in (3.34) and 

(3.35). As can be seen, the coefficient of e 3 involves only the Jj mode. 

2_ 

The vanishing of the coefficient of e 3 gives 


ib — 0 
v 0z 


(A. 27a) 


ai/) 2z - rp n*2 


/ Vie 3 C n - 1 2 ^ 

2 p il\ X lr^lr + -^ + X ll^l X iy 


(A. 27b) 


2 / 2 2 10 *10 

« X 2 z " rP ll X 2 = P 11 Kr - hr + ~2 T 


3 ^ll 1 2/2 2\\ 

+ 2 ^T1 ( X 1 “V 


(A. 27c) 
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The functions 0 q , ^ 2 > and X 2 are chosen to satisfy (A. 27a), (A. 27b), and (A. 27c). If 
i/jq is taken to be constant, (A. 27a) will be satisfied identically. Choose 0 2 and X 2 to ' 3e 


*2 ' E A On J o( X o„ r ) 
n = l 


c°sh[x 0n (z + h)j 

cosh X h 
On 




+ 7 1 A cos 2 8 + B sm29 J (x r) 

„ k 2n 2n / 2\ 2n / 

n = 1 


COSh [ X 2n (Z+h) ] 


2n 


(3. 37) 


and 


X 


-E 


2 i-i C On J o( X Oa r )' 

n = 1 


coshfx (z+h) 
t On 

cosh A h 
On 


w co sh[x 2 n (z+h) l 

/ (C cos 2 8 + D sin20)J (X r)— - 

i ' 2n 2n / 2\ 2n / cosh X h 


n = 1 


(3. 38) 


where 


J Vx a) = J '/x a) = 0. 
0 \ On / 2 \ 2n > 


By finding the appropriate generalized coordinates in 0 2 and X 2 , (A. 27b) and (A. 27c) 
can be satisfied. These generalized coordinates can be found by introducing (3. 34), 
(3.35), (3.37), and (3.38) into (A. 27b) and (A. 27c) and applying a Fourier-Bessel tech- 
nique using the following orthogonality conditions: 



0, m ^ n 



0, m 4 n 


2 2 

A a 

2n 


[ 2X 


2 

2n 



m 


n 


(3.39) 


(3. 40) 


These conditions give the generalized coordinates of the J 


0 


and 


modes as: 
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^0n 

’ °0nl 

( f l f 2 + f 3 f 4 ) 

A 

A 

2n 

°2n 

( f l f 2 - *3 f 4 ) 

A 

B 2n 

‘ n 2„( 

£ ! f 4 tf 2 f 3) 

A 

C 0n 

°0n / 

2 l 

f 2 - tf - *? - 

2 4 1 

Sn 

a 

"T$ 

.2 .2 , 2 
■2 +, 3 - f l - 

A 

°2n 

' n 2n( 

f 2 ' 4 - f l f 3 ) 


r (3.4i) 


j 


where O and Cl are constants defined in the appendix. 
On 2n 


The terms a r 


a ll b 00 


a ll b 01 


a 12 b 00 
+ + 


a ll b 00 b ll a 22 b 00 


00’ 01 « 02 a a ^ 
and each contribute to the coefficient of e. The coefficient of e contains sin cot, 

sin 2 cot, sin 3 cot, cos cot, cos 2 cot, and cos 3 cot. With this type of approximation it is 
assumed that only the first harmonic terms need vanish. The first harmonic terms 
from a^Q and are 


11 


rdx 


dr 


- v ip x 


r cos0 - r sin0 cos cot 


i i a * i \ 

J cos cot — + pX^sin cot 


/ 2 r 2r \ 

^ a - — cos 0 + — sin 0 j cos cot 


(3.42) 


2 2 
where p-^ r cos 0 cos cot corresponds to the translational motion uj, pp r sing cos cot 

corresponds to the translational motion u 2 , a cos cot corresponds to the translational 

2 r 2 r 

motion Ug, -^-cos 0 cos cot corresponds to the rotational motion c 02 > and — sin 0 coscot 

corresponds to the rotational motion co^. 


The first harmonic terms from a 0i 


a ll b 00 

are 

a 


11 


*lr X 2r - * MV 29 ' X l 9 *2 6 ) ‘ ' ’jfaVVSs) 


X ll«ll(*l X 2z - V 2 z) + ?(*1 X 2ZZ ' X 1*J 


cos cot 
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- p 


11 


X lr X 2r + ^lr^2r + 2( X 10 X 20 + ^ 1 Q^ 2 q) 


hl(hl - MW + *1*2) ~ X ll^ll( X l X 2z + *l*2z) 


7/x X„ +!M 0 ) Si: 

2\ 1 2zz Y \ Y 2zzJ 


sin cot 


(3.43) 


Similarly, there are first harmonic terms from a. 


02 


2a ll b 00 b ll + a 22 b 00 
2 or 2 


a ll b 01 a 12 b 00 
a 


which because of their length are not written out here. 


The equation obtained by setting the first harmonic terms of the coefficient of e equal to 
zero is now satisfied in a Rayleigh-Ritz, or averaged, sense by multiplying the equa- 
tion by 


r i( x u r ) 


cos 0 
sin 0 


rdrd0 


(3.44) 


integrating over the free surface, 0 < r < a, 0 < 0 < 2ir, and using the known results 
a 2 it 

1 cos 8 1 

r i( X li r ) drd9 = ° (3.45) 


0 0 


cos 0 
sin 0 


a 2 n 

I / 

0 0 

a 2ir 

If 


r sin 0 cos 


. 2 

sin 0 


9J i(\i r ) 


d r d0 = 0 


2 

cos 0 


0 0 

The contributions from aQQ + are 
a 2rr 


J i( x n r ) drde ' JL l J i( x ii a ) 

A-, 


11 


/ / ( a o« + B n) coseJ i( x n r ) rdrde = 


(3.46) 


(3.47) 


0 0 
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up 


11 


>/ df 2 

~ J l( x ll a )\lF 


2X 


V f. 


11 




11 


'i( x n a )] 


(«l) 


COS cot 


- up 


11 


. 2 2 
x ll a - 1 


2X 


11 


'i 2 ( x u a )(df + 


sin cot 


2u a 


'l( X U a ) 


(o) 2 ) 


cos cot 


(3.48) 


+ B n ^ sin 0J 1 ^X 11 r^rdrde 



= u p 


11 


v 2 2 

X li a il _1 2 

= J 


2X 


i( x n a ) 


11 


u a 


'll 


T l(\i a ) 


(u 2 ) 


cos cot - up 


11 


, 2 2 

X a - 1 „ 

— — — — — J 2 


2 X 


11 


i( x n a ) 


/ df \ 


+vi ) 

sin cot - 

\ dr 4 f 


2 u a 

X 2 h 
11 


’i( x n a ) 


(Wi) 


cos cot 


(3.49) 


The contributions from a 


01 


a ll b 00 


a 


B 2 are 


a 2u 



B 2 ooseJ i( x n r ) rdrde ‘ ’"ii J i(Vi) G i + f 4( f 2 f 


0 0 


V * 


/ 

V A , 

A 

)°2 

COS cot + fP 12 

f (f. 
2 ' j 

f j) G l - f 3 (< 

2 f 3 - f l f 4) G 2 


sin cot (3.50) 



f 4(Vih + h( £ 2 f 3 - hhK 


sin cot (3.51) 
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2 2 2 2 

where f. f. = f^ + f^ + f ^ + f^ , G^, and are constants defined in the appendix. 
The contributions from 


ab -ab 2abb + a b 
„ 11 01 12 00 11 00 11 22 00 

B = a - + — 

3 02 a „ 2 

2a 


are 


a 2 tt 


/ / B 3 o° s9Jl (x u r)rdrde = - ^ + f 4 ( f 2 f : 


0 0 


- f f ) 
1 4/ 


a 2 tt 


H 


cos cot - 


11 


A A 

f (f. f . \h - f (i f -f f )h 

2 \ j j / 1 3\ 2 3 14/2 


sin cot (3.52) 


j j B 3S ineJ 1 (x il r)rdrdS = - — ' ' f 2 ( f 2 f : 


0 0 


- f 


f 

1 4/ 2 


cos cot - 


11 


A A 

f 7f f \h + f/f f - f f \h 
4\ j ]/ 1 ll 2 3 14/2 


sin cot (3.53) 


where and are constants defined in the appendix. 


From the Rayleigh-Ritz process two ordinary differential equations are obtained. Set- 
ting the coefficients of sin cot and cos cot equal to zero in each of these equations results 
in four, first-order, nonlinear, ordinary differential equations. 

This system is 


df 

~d 


7= G i( f i-V VO- 1 = 1 - 2 - 3 ' 4 


where 


G 1 - “ H , 2 
G 2 = H ,1 
G 3=- H ,4 
G 4 = H ,3 


(3.54) 


(3.55) 
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and 


H = ( F l + A 2) f l + ( U 2 +A 2) £ 4 + ^Vi 

^^v/'Wvs-vJ 2 (3 - 56 > 

where 

H,i = ~~ (3.57) 

i 


The constants Fj, K , K , A , A^, U 2 are defined in the appendix by (A. 18), (A. 19), 
(A. 20), and (A. 21). 

A steady-state harmonic solution to the boundary value problem is given by the roots of 
the four equations 

G . (i , f , f , f ) = 0, i = 1, 2, 3, 4 (3.58) 

l \ 1 2 4 4 / 


The roots of (3. 58) are functions of fj where the f^ are independent of the time r. The 
form of (3. 54) is similar to the equations derived by Miles 7 for the undamped spherical 
pendulum. 

There are two solutions to (3.58). The first, called planar motion, is 


f l = * 


f 3 = 


f = f = 0 
2 4 


(3.59) 


where y is a parameter independent of time. The transformed frequency is 



- I^Ry 2 


(3.60) 
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where 


P - -(f 1 + a 2 ) 


Q = 


A 1 +U 2 
F 1 +A 2 


R = 1 + Q 

The second solution to (3.58), called nonplanar motion, is 


(3.61) 


f i =r 


2 V 7 K R 


f 3 = y<3 


1_ 
- 1\2 


Q - - LQ 


2 _ 2 Py" 1 

I — y 

4 7 K R 


= L 


(3.62) 


with 


-if A 2 K l\ 2 

" = -y K + — ) +R V 


2 / 


K 


K. 


K 


K 2 - 2K 1 


(3.63) 


(3. 64) 


It is seen that the nonplanar solution is real and, hence, exists for y > 0 when 


T 


KnR 


o ¥ 

> 0, and for y < 0 when y° - — — < 0. 

K2R 


The names planar and nonplanar motion are used in analogy with Miles' terminology for 
the spherical pendulum. It is not to be implied that the motion of the free surface is 
necessarily described by the names given to the two solutions of (3. 58). 
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SECTION 4 

STABILITY OF THE STEADY-STATE HARMONIC SOLUTION 


To determine the stability of the motion corresponding to a given steady-state solution, 
consider the perturbed solution 


f . . , (0) Xt 

f. (t) = f. + c. e 

1 i i 


c. « 1 

1 i‘ 


1, 2, 3, 4 


(4.1) 


The f^ 0 ) are constants corresponding to the steady- state amplitudes of the harmonic 
solutions of (3.54). The corresponding steady-state solution will be stable if Re(X) < 0 
and unstable if Re(X) > 0. 


Substitute (4. 1) into (3. 54), neglect products of the c^s, and use the fact that the f/ 0 ) 
are solutions of (3.58); the following set of homogeneous algebraic equations are 
obtained: 


d + X 

d 

d 

d 


c 


0 

11 

12 

13 

14 


1 



d 

d - X 

d 

d 


c 


0 

21 

22 

23 

24 


2 




d ™ 

d + X 

d 


c 


0 

31 

32 

33 

34 


3 




d „„ 

d 

d -X 


c 


0 

41 

42 

43 

44 


L 4 




where 


11 


2K f<°>f<°> + K f<°)f<°) 
1 1 2 2 3 4 


V + K (f. (0) f (0) ) + 2K f 
IV ] J / 1 : 


12 


( 0 ) 


- K f 
2 3 


( 0 ) 


d 13 - + - 2f2 <o, f3 (o, i 


(4.2) 


4-1 



j _ „„ , ( 0 ). (°) , ( 0 ) ( 0 ) 

d 14 " 2K 1 £ 2 U + K 2 f l f 3 




d = d 

22 11 


d = 2K f (0) f (0) + K f (0) f (0) 
23 113 224 


(0) f (0) of (0) f (0) 

24 “"1*1 *4 ' “21*2 3 


d_ , = 2K. f, ( 0 ) fj 0) + f^'f^' ~ 2f l ” ,f 4 


d = d 
31 24 


d = d 
32 14 


d = 2Kf<°>fW + Kf<V°> 

33 134 212 


d A = v + K 

34 1 


f <°> f <°> 1 +2 Kf <°> 2 - K f <°) 2 
j j 1 14 2 1 


d 41 d 23 


d = d 
42 13 


d = v + K (t ( 0 ) f (0) ) + 2K f 
43 IV j J / 1 


2 2 
<0) - v 2 (0) 


d = d 
44 33 


The solutions of (4.2) will be nontrivial only if the determinant of the coefficient matrix 
is zero. This condition gives an equation for the allowable values of 

4. 1 STABILITY OF PLANAR MOTION . Substituting (3. 59) into the expressions for 
the djj 's and expanding the determinant of the coefficient matrix in (4.2), one obtains 
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x 4 +X 2 (m i / + M ' + 2M 3 J + M^Mg + 4K x 2 y 4 Q^M 4 - K 2 Q 2 y 4 M 5 + 


M„ 


where 


M ' 


v + y 2 (n 1 R - k 2 q) 


M 2 


p+K^y^l+3Q^ 


= 


K 2 y Q 


M. = 


M = 


v + y 2 (K : R - K 2 q) 

2 / 2 \ 

V + K y ( 3 + Q J 


V + K^y ^3 + Q J 


v + y 2 ( k 2 r - k 2 ) 


2 y Q 


" + A K 1 R - \) 


2 / 2 \ 
V + K y (l + 3Q ] 


The boundary between stable and unstable planar motion corresponds to X = 0. 
X = 0 in (4.3) and substitute forp from (3.60): 


jN 4 y 12 - Py 9 + N^y 6 + P^y 3 + P 4 


= 0 


= 0 (4.3) 


Set 


(4.4) 
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where 


N 1 ‘ - K 4 + 2K 1 Q2 - K 2 Q 

N 2 = K 2 K 4 Q - K 2 Q 2 (K 2 * 2Kj - 2 Kl K 2 (l + Q 3 ) 

N g = 2K 1 K 2 2 (q 4 - Q 3 + Q 2 - q) 

2 2 2/ 2 \ 

N 4 = -rK 2 K x Q (Q + 1 ) 


(4.5) 


One possible solution to (4.4) is y = ±oo. Since y is actually an amplitude this would 
correspond to unstable motion. Letcr = y 3 ; then the other possible solution to (4.4) is 


4 a o 2 3 4 

No - N Pa + N P a + P N a + P -0 
4 3 2 . l 


(4.5) 


This equation can then be solved for a = y 3 . The solution of (4. 5) for various pertur- 
bations is given in Section 5. 

4.2 STABILITY OF NONPLANAR MOTION . Substitute (3. 62) in the expressions 
for the dy 's: 

d n = yQK 4 L 


K Py 


-1 


d 12 = ^ + ’' 2 ( 2K 1 R - - Tr(‘ + 3Q2 ) 


13 


- Ly(K 4 Q 2 + K 2 R) 


/ 2 2K 1 P ^ 
d^ = Q(K 4 y + 


- 1 . 


K 2 R 


21 


v + y 2 ( 2K i R ■ k 4 ) + p r -1 ^ 




K 


d 22 ’ d ll 


23 


Q Ik 4 y 


2 P- 


R 
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r 


d„„ = -yL^K, + K R) 
24 r \ 4 2 / 


d = d 
31 24 


d„„ = d., _ 
32 14 


33 


= -d 


11 


-1 


d 34 ‘ «-r 2 (2K lR -K 4 ) 


K Py 

-hr(-« 2 ) 

2 


41 


= d 


23 


d 42 d 13 


d 43 = V + 


: ( 2K i R - K 4 « 2 ) - P >'" 1 (- 


- 1/ K 1 Q 2 ' 


d = d 
44 33 


Substitute the value of v given by (3.63) into the expressions for d_^> d o ^, d , d : 

2 -1 

2K PQ y 


2 / \ 1 

d = y (2K R + K ) - ; 

12 r \ 1 4/ ] 


K R 
2 


d 21 = r 2 (2K lR + K 4 Q 2 ) + 


-1 


R 


34 


2 K Py 

! (2K lR + K 4 Q 2 )- 


R 


d 43 ” y ( 2K i 


* + K 4 ) + ^ 


Substitute the above values of the dy 's into (4. 2) and set the determinant of the coeffi- 
cient matrix equal to zero to give a fourth-degree polynomial for the determination 
of the parameter X, such that (4.2) has nontrivial solutions. Stability of the 
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steady-state nonplanar solution is determined by examining the roots of this quartic 
equation. The regions of stable and unstable motion are given in Section 5 for various 
perturbations . 
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SECTION 5 


NUMERICAL EXAMPLE 


To compare results obtained here with those obtained by Hutton, use the following values 
of the parameters: 


a = 

tank radius = 5. 938 inches 




h = 

water depth = 8.907 inches 




X 11 

a = 1.84119 


► 

(5.1) 

J l(x n a) = 0.581865 

> 



Then, 





F 1 

= 8.53992 




e 

= 0.99205 




p n 

= 10. 987 rad/sec = 1. 734 cps 




A 2 

= -0.6389 


► 

(5.2) 

A 1 

= - A 2 




U 2 

= F 1 

- 



To evaluate and K 2 , use only the first five terms in the 

and G 2 to approximate the series. In the calculation of 

infinite series defined by 
Gi and Go the last three 

terms are about one percent of the first two terms. Thus, 




K 1 

= 0.4853 X 10 5 




K 2 

-4 

= 0.13707 X 10 

! 

► 

(5.3) 

K 3 

= -3.0235 




K 4 

-6 

= 4.0010 X 10 
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5.1 PLANAR MOTION . Equation (3.60) gives the transformed frequency, p, in terms 
of the parameter y. The coefficients in this equation depend on the perturbations given 
to the liquid-tank system. Thus, 


V 


- ( F i + \Y' 


- K. 


+ 




(5.4) 


Equation (4.4) gives the values of y which separate stable and unstable regions. The 
coefficients in this equation also depend on the perturbations given to the liquid-tank 
system. 


5.1.1 Case 1. Consider ajj = to 2 ~ u 2 = 0, Uj = e cos cut, which is the case con- 
sidered by Hutton. Then, 

-1 -5 2 

V = - 8.5399y - 0.4853 X 10 y (5.5) 

The motion is unstable for 


- 0.1337 < p < 0.06459 


(5.6) 


where 


y = 95.82 when v = -0.1337 

and 

y = -85.41 when p = 0.06439 

3 

This case agrees with the result obtained by Hutton . 

5.1.2 Case 2 


(5.7) 

(5.8) 


C0„ = T" COS COt 

2 h 

-1 -5 2 

V = 0.6389y - 0.4853 X 10 y 

The motion is unstable for 

- 0. 0237 < v < 0. 0115 
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where 


y = -40.377 when p = -0.0237 


and 

y = 35.988 when p = 0.0115 
5.1.3 Case 3 


u 


1 


0 


- <X > 2 = — cos cot 


“1 “6 2 
V = 0. 6389 y - 0.9706 X 10 y 


The motion is unstable for 


- 0. 1863 < v < 0. 1934 


where 

y = -3.43 when p = -0.1863 

and 

y = 3.30 when p = 0.1934 
5.1.4 Case 4 

e 

to. = — COS cot 
i h 

W 2 = u 2 " 0 
= f cos cot 

-1 - 5 2 

1/ = - 8. 5399 y - 0.488 X 10 y 

The motion is unstable for 

- 5.735 < V < 5.662 




> (5.9) 


(5.10) 


A 


(5.11) 


(5.12) 
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where 


y = -1.508 when V = -5.735 

and 

V = 1.489 when j/ = 5.662 

5.1.5 Case 5 

to = u =0 
1 2 

e 

a) = — COS COt 
2 h 

- e cos cot 

V = -7.901 y 1 - 0.4853 X 10~ 5 y 2 
The motion is unstable for 

- 0. 1269 < V < 0. 06133 

where 

y = 93.37 when v = -0.1269 

and 

V = -83.23 when v = 0.06133 

5.1.6 Case 6 

= u = e cos cot 

^ = w = -| cos cot 

-1 -5 2 

V = - 7. 9010 y - 1.144 X 10 y 

The motion is unstable for 

-4.772 < v < 4.834 


(5.13) 


(5.14) 


(5.15) 


(5.16) 
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where 

y = 1.655 when v - -4.773 

and 

y = -1.634 when u = 4.834 
5. 1. 7 Case 7 


m = ~ cos wt 
1 h 

u = e cos oot 

Li 

The results are the same as in Case 5 (by symmetry). 
5.1.8 Case 8 


to. = 4- cos cot 
l \ h 

u„ = e cos cot 
2 

The results are the same as Case 4 (by symmetry). 

Cases 1 and 2, according to (2), should not differ. However, it is seen that the unstable 
region for Case 2 is much smaller than the unstable region for Case 1, indicating that, 
at least for stability considerations, rotational oscillations about x 2 are not equivalent 
to translational oscillations in the direction. Also note from Case 5 that the unstable 
region is slightly smaller but nearly equivalent to the unstable region in Case 1. Here 
the rotational and translational motions are taking place in the same plane. The rota- 
tional motion thus has a much smaller effect on the free-surface motion than does the 
translational motion, even to the extent that the combination is essentially not different 
than the situation for translation alone. 

Cases 4 and 8 consider the combination of rotational motion and translational motion in 
planes perpendicular to one another. From (5.12) the region of unstable motion in 
these cases is much greater than any of the other cases considered. 
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5.2 NONPLANAR MOTION 


5.2. 1 Case 1 


CO, = co„ = u„ = 0 
12 2 


u = e cos cot 


The quartie equation for this case is 


x + (m 3 + mJa + m 5 m 6 = 0 


where 




d d 

34 13 


d d 

24 12 


M = 


12 


24 


13 


21 


M = 
5 


d d 

13 12 


d 12 d 13 


M = 


24 


34 


21 


24 


and the d 's are defined in (4.6) and (4. 7). 
The transformed frequency, u, is 

/ A "I 

v ' -( F i + y +E V 

2 

For this case, 


A 2 = ° 
R = 1 


F = 8.53992 


_ 1 o 

v = -3. 0235 y + 4.001 X 10 y 


(5.17) 


(5.18) 


(5.19) 


(5.20) 


(5.21) 
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The generalized coordinates f, for this case are 

(5.22) 

> 

Since must be real, f > 0. Thus, y cannot be in the range 

F 

(-85.41) = - < y 3 < 0 (5.23) 

Evaluate M + M and M M : 

3 4 5 6 

= 7.5148 x 10 y (y* + 3.5699 x 10 5 j (5.24) 

M 5 M 6 = 2 - 5677 * 10 15 y 1 (y 3 + 6.2305 X 10 5 ) (y 3 + 3. 7784 X 10 5 ) (5.25) 

2 

Consider (5.17) as a quadratic in A. : 

2 

(A 2 ) + (M 3 + M 4 )( X 2 ) + M s M g = 0 (5.26) 

If 

3 5 

y <-6.2305 X 10 

or 



y < - 85. 41 

M + M, > 0 
3 4 

M M < 0 
5 6 

Thus, by Descartes' rule of signs there is one positive real root of (5.26). This corre- 
sponds to a region of unstable motion. The case - 6.2305 x 10 5 < y 3 < 0 is dealt with 
in (5.23). 
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To examine stability for y > 0, compute the discriminant of (5.26): 

A = 52.99 X 10 -2 °y 9 + 361.4 X 10 14 y 6 + 572.5 x 10~ 10 y 3 

- 241.75 X 10” 5 (5.27) 

20 3 

Set A = 0, multiply by 10 , and let a = y to obtain 

52.99cr 3 + 36.14 X 10 5 a* + 572.5 X 10 10 a - 241.75 X 10 15 = 0 (5.28) 

The only positive real root of (5.28) corresponds to y = 64.47. Thus for 0 < y < 64.47, 
the roots of (5.26) are complex. If y > 0, Mg + and Mg Mg are both positive. Thus 
by Descartes' rule of signs there is no positive real root of Equation 5.17. Replace 
y 2 by -a in (5. 26) to obtain 

01 ~ ( M 3 + M 4)°' + = ° (5.29) 


Thus by Descartes' rule of sign there exist either two positive real roots or no positive 
real roots of (5.29). If y > 64.47, the roots of (5.29) are real. Hence for y > 64.47, 
Re (y) ^ 0, and this relation corresponds to stable motion. 

In summary, the steady-state nonplanar motion is stable for 


64. 47 < y < oo 
- 0. 03027 < p < 00 
and is unstable for 

0 < y < 64.47 
-co < i) < -0.03027 
- co < y < - 85. 41 
0. 06459 < p < co 


(5.30) 


(5.31) 


(5.32) 


The solution does not exist when 


-85.41 < y < 0 

-co < p < 0. 06459 
2 

since then f^ <0. 
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This motion is stable for a small range of driving frequencies, co, which includes the 
first natural frequency, Pjj. This condition can be seen from the range of transformed 
frequencies, u, given by (5.30). 

O 

The results have agreed with those obtained by Hutton 0 . 

5.2.2 Case 2 


go. = 4- cos oot 
2 h 

Proceeding as in Case 1, with = 0 and = - 0.6389, the quartic equa- 

tion to be examined is 

4 / \ 2 

X + (M + M X + M M = 0 (5. 33) 

\ 3 4/ 5 6 v ' 


where M3, M4, Mg, and Mg are defined by (5.18) and (5.19). Here the d-y 's have dif- 
ferent values than in Case 1. 


The transformed frequency, v , is 
K 


_1 -1 
^2 


- a 2TT’' + V 


or 


— 1 —6 2 
= 0.2262 y + 4.001 X 10 y 


(5.34) 


(5.35) 


The generalized coordinates for this case are 


f i = y 


f 3 = 0 


2 _ 2 2 -i 

f 4’ y + V 


(5.36) 


Since f must be real, > 0. Thus, y cannot lie in the range 


A 3 2 

0 < y < - ¥ - 
2 


(5.37) 
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where < 0. This restriction gives 
0 < y 35. 8 
0. 01047 < v < co 

2 

Regard (5. 33) as a quadratic in A and examine the roots of 


(5.38) 


(a 2 ) 2 + (m 3 + mJ a 2 + m.m 


5 6 


0 


(5.39) 


This gives the following regions of unstable and stable motion, respectively: 

-co < y<0 and -» < v < « (5.40) 

35.8 < y < co and 0.01047 < v < 03 (5.41) 

This motion is not stable about the first natural frequency, as opposed to the situation 
in Case 1. 

5.2.3 Case 3 

w, = 0 


60 „ = 


u„ = 


— cos cot 
h 


e cos cot 
0 


The quartic equation for this case is (5. 17) with M 3 , M 4 , M 5 , and Mg given by (5. 18) 
and (5.19). The d^ 's are evaluated from (4.6) and (4.7), with 


A, = U = 0 
1 2 


Q = 0 
R = 1 


The transformed frequency, p, is 

“1 —6 2 
V = -2.797y +4. 001 X 10 y 


(5.42) 
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The generalized coordinates are 


f i = y 


f = f = 0 
2 3 


f 2 _ 2 . 
f 4 - y + 


( F i +A z) -1 


K 


Since is real, f > 0. Thus, y cannot lie in the range 


/F, + A \ 

< y 3 < 0 


(5.43) 


(5.44) 


Thus, the solution does not exist for 
-75.6 < y < 0 
-o= < v < 0.0597 


(5.45) 


The regions of unstable motion are 

-co < y< —8. 31 and 0.0613 < V < °° (5.46) 

0 < y < 62.58 and - m < i> < - 0.0399 (5.47) 

The region of stable motion is 

62.58 < y < oo and - 0. 0399 < v < co (5.48) 

Note the similarity between Cases 1 and 3 with regard to the stability of the nonplanar 
motion. The rotational motion has a much smaller effect on the stability of the free- 
surface motion than does the translational motion, even to the extent that the combina- 
tion is essentially not different than the situation for translation alone. 

In all cases considered for the stability of nonplanar motion, there is a region of p for 
which it is possible to have both stable planar motion and stable nonplanar motion. 

5.2.4 Case 4 


u 


2 




0 
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U) = 4- COS COt 

1 h 

= c cos cot 

The quartic equation for the determination of the allowable values of X is 
4 2 

X + MX + N = 0 

where 



d 

d 


d 

d 


d , 

d 


11 

12 

+ 

11 

12 

+ 

11 

12 


d 

d 


d 

d 


d 

d 


21 

22 


32 

33 


43 

44 


d 

d 


d 

d 


d 

d 


22 

23 


22 

23 


33 

34 

+ 

d 

d 

+ 

d 

d 

+ 

d 

d„„ 


32 

33 


43 

44 


43 

44 


and N is the determinant of the coefficient matrix in (4.2) with X = 0. 

The transformed frequency, v>, is 

—1 —6 2 

V = -3. 0235 y + 4.024 X 10 y 

and the generalized coordinates are 

(2 5 -1] 2 

f = - 0. 0748 | y + 6.24 X 10 y J 

j 

•> 

f = 0. 0748 y 

3 ' 

2 2 5 -1 

f = y + 6.24 X 10 y 

4 

2 

Since f is real, f >0. Thus, y cannot lie in the range 
4 4 

-85.41 < y < 0 and 0.06459 < v < <*> 


(5.49) 


(5. 50) 


(5.51) 


(5. 52) 


(5. 53) 


(5. 54) 
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The motion is unstable for 


-ra < y < -85.41 and 0.06459 < t/ < <*> 

59. 1 < y < co and - 0. 0371 < v < 03 
The motion is stable for 

31.6 < y < 59.1 and -0.0915 < V < -0.0371 
This case, as in the planar motion, has the largest region of unstable motion. 
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(5.55) 
(5. 56) 

(5. 57) 


I 


SECTION 6 


CONCLUSION 


This paper considers the irrotational motion of an incompressible, inviscid fluid con- 
tained in a partially filled tank. The tank is subjected to both transverse and rotational 
vibrations whose frequencies are near the first natural frequency of small free-surface 
oscillations. The analysis was performed using a method originally suggested by 
Hutton^, retaining higher-order terms in the free-surface dynamic and kinematic bound- 
ary conditions. The theoretical investigation predicts the forcing frequency ranges, for 
various combinations of rotational and translational motion, over which there are stable, 
steady-state, harmonic, planar and nonplanar motions. The least stable case occurs 
when a combination of motions occurs in planes perpendicular to one another. This 
condition substantiates the findings of Hutton in that the mechanism that apparently 
causes sloshing in the unstable regions is a nonlinear coupling of the fluid motions par- 
allel with and perpendicular to the plane in which the translational motion is taking 
place. 
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APPENDIX A 


In the sloshing problem being considered, the free-surface height, t?, is an unknown 
which may be eliminated by replacing (3. 24) and (3. 25) by a single equation which does 
not contain tj. Solving (3.24) for?}, one obtains 

av = r(r, 0, t; (r, 0, t), t) (A.l) 


where 


“ T^r, 0, r] (r, 0, t), tj = + [<* - ra^ sin 0 + ra> 2 cos 0 - a^u.^ 


+ tOi u 2 + oo 3 oo 2 r sin 0 + co 3 r cos 0 


/("] 1 + "2 ) ' 


a - rci sin© 


1 2 

+ ru 2 cos0 - cOgU^ + OJ 1 u 2 + co 3 u)^ r sin0 + to 3 co 1 r cos 0 j 

2 

( 2 2 \ 1 ( 2 *^0 2 \ 

W 1 + W 2 j i\f r + + K) “ h + V COsS + a 2 rsin9 +U 3 3l 

u 2 

32/22 2 2 \ / 3\/ \ u 3 

+ 2 r (co 1 sin 0 + 00 2 cos QJ + 2 ^ + —) (ru^ sin0 - rto 2 cos0j + — 


+ w 3 r ^ sin0 - u^ cos 0 j - 3 r 2 oo 2 sin0 cos 0 + <x> 3 1 /)^ 


K + 


and the negative sign in the quadratic formula is chosen so that ari remains finite if CO ] , 
ou 2 , Wj , co 2 are all equal to zero. Equation (A. 1) can now be used to obtain the partial 
derivatives, rj t , r) r , and rjg. Thus, 


(« - r „h - r t 

(“ - r >r * r r 

(“ - r „H ' r e 


(A. 2) 


A-l 



Multiply (3.25) by (a - F^) and use (A. 2): 

(~4> z " u 3 + 2r to 2 cos6 - 2 rto i sinej ^q; - r ' 'j 

= rt-rA-^V^sjV. “ n z = ” - (A - 3> 

The potential functions in (A. 3) are evaluated on z = 77 , and, thus, (A. 3) depends upon 
T) implicitly. However, (3.24) depends upon 77 both explicitly and implicitly. The wave 
height, 77 , is eliminated between these two equations by first expanding the functions 
defined by these two equations in Taylor series about z = 0. 

Introduce the notation 


m n p s 
9 z v t 



m n p ,s 

ar ae az at 


z 


0 


(A. 4) 


where 


m + n + p + s = k 

The Taylor series of the function defined by (A. 3) and (3.24) can then be written in the 
form 


a 2 2 a 3 3 

a 0 + a j*7 + JT 7 ? + JTV + = 0 


(A. 5) 


and 


b 2 2 b 3 3 

b o + V + iT 77 + ¥ ! 77 


+ . . . =0 


(A. 6 ) 


2 . 2 


respectively, where 

b 0 " - b oo - V - V sin9 ■ + “2 si ”^) 

2 

U„ „ Ur. 

/ 3 \ 2 . o 

- 2 £ ^rco 1 sin 0 - r 00 ^ cos 9 + -j" j + 3 r u^ 00 2 sin 9 cos 9 + — 


- u rw sin 9 + u uj r cos 0 - u co r sin 9 + co u r cos 9 - to £ 

O 1 o Z io o Z Or? 
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b ± = -b u - b 12 - a + rw 1 sin0 - o^rcose - 


sinQ 


3 

- ra^cose + yj + W 2 u 1 - oj 3 oj 2 r sinQ -to 3 4 


8z 


- a3 3 co i rc °s8 - 


b 2 = - b 22 - b 23 + (“l + " 2 2 ) ' 2 ^z( r “l S “ 8 ' ra, 2 COS ® + t) 


- 03 £ 

3 ’0zz 


b„ = -£ + u, r cos 8 
00 s t 1 


, 1 I 2 ^0 2 

01 2 Pr 2 5 z 


3 b 


b, „ - 


00 


11 3z 


3b 


01 


12 gz 
3 b 


00 


00 St 


3b 


mn 


m+1 , n+1 


3 z 


for m, n = 1,2,3... 


(A. 7) 


It is evident that the potential functions are of the same order as the wave height, as 
can be seen by neglecting the products of rj and £ in (A. 6); then, the first approximation 
becomes 

-b 00 -otn = u 2 r sin 8 


or 


V 




r cos 0 - u r sin 0 


(A. 8) 
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With this fact in mind, one can expand the function defined by (A. 1) in a binomial expan- 
sion and neglect terms of 0 (rft)- Here it is assumed that ^3’ u l> u 2’ u 3’ ant ^ 

their time derivatives are of the same order as the wave height. This assumption gives: 


r(r,0,T)»t) = 


1/2 21 

r cos 0 + — + — + / + ^2 r sin ® + u 3^5 


3 r 2 / 2 2 2 2 \ / u 3 \ / \ u 3 

+ — — (cOi Sin 0 + C0 2 COS 0 J + 2^ + — J (ru^ sin0 - rw 2 cos QJ + — 

+ co^ r^sine - ^cosflj - 3 r 2 sin© cos0 + 

+ -I- ro) sin0 + roi)„ cose] 1-0 + u r cos0 + u r sin0 
all 2 Jit J. £ 

( 2 2 
2 ^0 2 j 3 r 2 / 2 .2 2 2 \ . U 3 

+ ~2 + ^z) + TK Sm 9 + a, 2 C ° 8 9 / + U 3^ + T 


+ 2 sin0 - rco 2 cosdj + u> 3 r ^ sin0 - u 2 cos^j 

- 3 r 2 C0 1 U) 2 sin 0 cos 0 + + 0 1 1” ^2 U 1 + W i U 2 + u5 3 r (w 2 sin0 

\lf • • . „1 -2 f 2 . 2 .2 

+ a) cos 0 ) + r cos 0 + u 2 r sin0 - a [ r w i sin ® 


+ r 2 a; 2 cos 2 0 - 2 r 2 uj^ u:^ sin 0 cos 0 j | - 0^ + u^r cos 0 + u^ r sin 0 j 

+ 0(7? 4 ) ( A - 9 ) 

Use (A. 9) to find r , r„, and T . Note that on z = rj, T = T \ then o- - T can be 
rot 77 z Tj 

computed. Substitute these values into (A. 3) and expand the function defined by the 
result in a Taylor series about z = 0. This operation leads to (A. 5), where a^ is 
defined below: 

for k = 0, 1 , 2, 3 . . . 

z = 0 
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where 


© = 


*e*et 


- a(i z - </>„ + % r oosS + 2* r (S rt + 2 — j- + 2« z * zt * “ 3 “3 


- ^ (♦j.oobS -tg 4^* 


*8*86 ,2, _ 

r rr ' 4 " *z*zz + 3 


*r*8 


- 2lb lb lb - 2 

Y r Y rz Y z 2 


*r*8*r8 „ *z*8*8z ( , | 

— + (- U 3 + 2r | 


- 2 


0J 2 COS0 


1\I Ve« / 

- Wl smej) a + >p r ^rz + — — + Mzz ' *tz ' 2r ^zz f 2 COS0 


r^rz 2 

r 


1 smej + u 3 !() zz tu 3 l(i ez - a' 1 r(i 2 cose - « 1 »*»8)/* r # rz 


, *0*0 z , , 

-*tz + ^- + *z*zz 
r 


-^ z [- 2 ^zz( U, 2 COS0 - W l Siae ) + W 3^e z 


Mrt 


+ u 3 i/) zz +a 1 r^to 2 cos0 - 0)2 sine^i/j^ j -a 1 r^a) 2 cos0 -<I) 1 sin0^ 

i/j il) 

# 0/ .2 2 

+ + 0 ih - 0 +u r cos 0 + u„ r sin0 + 3 r* 1 [cu.ci) sin 0 + a}r,k>„ cos 0 

2 z r zt tt 1 2 \ 1 1 z 2 

r 

- 2^[) zt + -^-^r^co 2 cos0 -a) 1 sin©^ - 2r ($> z + cos ® ^ sin0^ 

+ sin0 -U 2 cos0^ + co 3 r ^ sin0 - u 2 cos0^ + u 3 u 3 + u 3 i/> zt + u 3 i/) z 

- 3 r 2 sin0 cos 0 a) 2 + w 1 a) 2 j + u> 3 ip^ + oo 3 ipg t 


u r sin0 + u„ ip 

2 <j 2 


+ sr 2 ^^ sin 2 0 + oo 2 ob 2 cos 2 0j - 2^) zt + y^r |u> 2 cos0 -co x sin0^ 
- 2r + y^ (ui^cosd - ty sin0^ (u.^ sin0 - u 2 cos0j + Ugi/J 


A -5 



i / , 2 % , 2 ) 


-0 t + i^cose + 2\^ r + “ + 0 , 


+ Wgr^u^sinQ - i^cosej - 3 r 2 sin0 cos 0^0^ C0 2 + u^ci^j + ^3^1 
+ C0 3 *^0 t ] - a cosQ - c^ sin dj 

+ r sin 0 + -^r^co^ sin 2 0 + co 2 cos 2 ej - 2r|^ + -j- ) ^ cos 0 

2 

\ / \ 2 U 3 

- co. sin 9 ) + w„ r(u sin0 - u cos 0) - 3 r to, m sin 0 cos 6 + u 0 tb + 

1 / 3 \ 1 2 / 12 3^z2 

- + Wl u 2 + V (u 2 slne * WjCosS)] [- fe 


+ 00 ll) 
3*0 


+ r cos 0 + r sin 0 - a 1 _ ^ 2 U 1 _ ^2 ^1 + ^1 U 2 + 00 1 S 


- 2r 


CO 3 r^to 2 Sin0 + co^ cos 6 j + CO3 r^co 2 sin0 + cos0j j [ ~$ t + ^ r cos0 

1 -2f 2 . .. 2 2 . .. 2 

u 2 rsin0 + « I 2 r sin 0 + 2 r co 2 co 2 cos 0 

sinecose^ai^a^ + ^ 2 )] [ + ^1 rcos ® + r sin ®j 

-2( 2.2 2 2.2 2 2.. if 

+ a I r co sin 0 + r co cos 0 - 2 r CO co cos0 sin© -0 +li r cos0 
L-L Z J. Z J L tt I 

if. / 2 2 2 2 v / 

+ u r sin0 I + ij) -u sin0 - 3 r(co sin 0 + co cos 0 ) + 2ip r cos 0 
z j r z \ j. z j z i? \ z 

- a^sin©^ + 2 ^ z + "^)( t °2 cos 0 " co i sin©j - sin© - u 2 cos 0 j - U 3 0 

+ 6rco i co 2 sin0cos0 - + a ^ci^cose - d^sinej (~ip Tt 

+ cos0 + u 2 sin0 j + a "^co 2 cos0 - co^ sinej + rcos© 


u 2 r sin© j j + co 3 |j/>g t +u ;L rsin0 - i^rcose - 


Wee 

Mflr “T" - U 3^0z 
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2 2 2 2 / 

- 0 0 - 3 r CO- sin Q cos0 + 3 r uj cos0 sin0 +20 r(c 0 o cos 0 

2 0Z X 2 Zq \ 2 

- u^sinfij - 2r ( l / ) z + '^)( a5 2 Sin ® + k^ 0015 ®) ” W 3 r^u l cos0 + u 2 sin0^ 
+ 3 r 2 a) 1 co 2 ^cos 2 0 - sin 2 ©j - co 3 0 ge + a 1 r^co 2 cos0 - ci^sinoj ^-0 

- u^rsinO + u 2 rcos0^ - a 1 r^co 2 sin0 + to^ cosO^-^^ + u^rcos© 


1 u r 2 2 2 2 

+ u 2 rsin0j + — - u 2 r cos 0 - 3r u^sinOcosS + 3 r co 2 cos0sin0 

r U 3 

+ 2 ^ z q r ( c0 2 c°s0 ~ aJ^sinej - 2r (^ z + 2 _ )( aJ 2 Sin ® + cos ®) 

- to r(u cos6 + u sm0)+ 3 r 2 a) to (cos 2 0 - sin 2 0) - t 0 o 0 Q . - u 

O'-L 2 / 1 2 \ / 3 00 «j 2 

+ a r (co 2 cos© - to^ sin0 j ^-0 " r sinQ + u 2 rcosQj 

- a r (^ 2 sin ® + cos ®) ( + a l r cos ® + ^2 r sin ®) J + °(^ 4 ) 


Equation (A. 6 ) is solved for 77 : 


b / b \ 

0 / 2 \ 2 

” = ^ w + 


0 n/ 3 \ 
= + 0(77 ) 


(A. 10) 


since b 2 = 0 ( 77 ) . Substitute (A. 10) into (A. 5): 


a + a I - 
0 1 


b a b 2 

i[-b7 + °( T » )J + ~ [“b^ + °( T » )j + ) 


a + a 
0 1 


a /b 2 
2 / 0 


(A. 11) 
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Compute the indicated multiplication; (A. 11) is then 


where 


ab ab ab abb 

11 00 , 11 01 12 00 11 11 00 

a + a _ — + a - + + 

00 01 q; 02 (x o; 2 

a 

a b 2 , 

22 00 „ „ / -1 . . „ -2 . . 

+ _ + B n + B l2 + B l3 - a n / a u 2 rsmQ + a raising 

2 ! a \ 

-2 . „ -2 . — 1 f 3 2 2 2„ 3 2 2 2 

- a b^ u g r smfl + a rcr 1 b 00 + a 2 ra5 i sm9+ 2 r COS 0 

+ 2 ^r( U 3 “ a i) " 3 w 1 a) 2 sin 0 cos 0 + Su^c^ + ^ cursing +co 3 £ Q \ 

- ^(“'>00 + V 8 *” 9 ) - - roj[b 00 + i, 2 rsm9] 

-if, 3 p 2 2 3 o 1 2 / v 

+ “ [ b 01 + i r “l sm 9 + i r *1°°* 9 + 2 S 2 r ( U 3 - °l) 

“3 9 e]) + a l2 a ' 1;i : 


- 3r u) co sing cos 0 + 2u ra + u CO rsinQ + 
12 3 113 


r sing 


Q 

-If 1 22 -2 / 2 2 2 \ 

+ B 23 a [ b 00 +i 2 rsin9 J + — “ ( 2b 00V Sine + \ r Sbl V 


33 -2 [ . 12 / 4\ 

T7 a [ b 0 o + u 2 rsine l + o(v ) = o 


(A. 12) 


00 


01 


02 


' a «z - *tt + V COS ® 




2£ £ + - 9 9t + 2£ £ - u/£ cos 0 - £.-^.\ 

rt 2 s z s tz l\ s r s 0 r / 


sin0 




1 2 t 
’0 *00 


r rr 


£ £ 

2 9 r *0 

- £ £ + — -2- -2i f £ 

z ^zz 3 r ^rz z 


2 ^r^0^r0 2 ^z^0^0z 
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a il -a ^zz " ^ttz 


3a 


mn 


m+l,n+l 3z 


ct = cos0 - sin0 

A- Ci 1 


_ 3(T 1 

r i,e Be 

Scr , 


CT, = 


i at 


j8 = u sin0 - u cos0 

-LI Ct 


P 


1,0 30 


0, = 


1 Bt 
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In summary, the boundary value problem in terms of the potential function, with the 
higher order approximation for the free-surface condition, becomes 
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and (A. 12). 
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Here p is the lowest natural frequency of small, free-surface oscillations; 
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